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ABSTRACT 

The  prediction  of  turbulent  secondary  flows  with  Reynolds  stress  models  in  circular 
pipes  and  non-circular  ducts  is  reviewed.  Turbulence-driven  secondary  flows  in  straight 
non-circular  ducts  are  considered  along  with  turbulent  secondary  flows  in  pipes  and  ducts 
that  arise  from  curvature  or  a  system  rotation.  The  physical  mechanisms  that  generate  these 
different  kinds  of  secondary  flows  are  outlined  and  the  level  of  tiirbulence  closure  required 
to  properly  compute  each  type  is  discussed  in  detail.  Illustrative  computations  of  a  variety 
of  different  secondary  flows  obtained  from  two-equation  turbulence  models  and  second-order 
closures  are  provided  to  amplify  these  points. 
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1.  INTRODUCTION 


The  turbulence  structure  of  internal  flows  within  circular  pipes  or  non-circular  ducts  can 
be  altered  considerably  by  the  occurrence  of  secondary  flows  [1-21].  Tlu'se  sec  ondary  flows 
lead  to  friction  losses  and  can  shift  the  location  of  ina.xiinum  inoim'iituni  transport  from  the 
pipe  or  duct  centerline  two  effects  that  can  have  profound  consetpiencc's  for  ('iigineerinp; 
design.  C’onsequently,  there  is  the  need  for  turbulence  models  that  can  ix'liably  |)redi(  t 
the  secondary  flows  that  occur  in  engineering  applications  which  include'  turbomachinery 
impellers  and  blade  passages,  aircraft  intakes,  and  |)ii)e  or  duct  cooling  systems,  to  name'  a 
few. 

In  this  paper,  the  prediction  of  four  fundamental  tyjies  of  secondary  flows  are  discusse'd: 
(1)  turbulence-driven  secondary  flows  in  straight  duets  e)f  neni-circidar  ereess-se-et ie)n,  (2) 
turbulent  se'condary  flows  in  curveel  circular  pipe's,  (3)  turbuh'nt  se'eeenelary  ne)ws  in  curve'el 
ducts  of  non-circular  cross-section  and  (d)  turbulent  .st'conelary  fle)ws  in  re)tating  ducts  eef  ne)n- 
circular  cross-section.  These  flows  are  se'lecte'd  sinee  thew  involve  se'ce)nelary  flows  ge'iie'iate'd 
by  a  combination  of  the  effects  of  normal  Reynolds  stress  elifference's,  streamline'  curvature' 
and  body  forces  arising  from  a  system  rotatieui.  Thus,  a  re'lative'ly  bre)ael  basis  fe)r  the' 
evaluation  of  models  can  be  provich'd.  The  ability  of  twe)-eeiuatie)n  mo'lels  anel  se'ce)nd-e)rele'r 
closures  to  predict  these  type's  of  turbulent  se'conelary  flenvs  will  be  e'valuate'el  in  a  systc'matie- 
manner.  A  variety  of  illustrative  calculatieuis  of  secondary  flenvs  will  be  pre'senlenl  along  with 
an  assessment  of  the  progress  that  has  been  made  in  the  analysis  e)f  tlu'se  flenvs. 


2.  ANALYSIS  OF  SECONDARY  FLOWS 

We  will  consider  the'  mean  turbulent  fleiw  of  a  vise  ous  inceunpre'ssible'  fluiel.  d'lie'  Re'vnejlels- 
averaged  Navier-Stoke's  and  continuity  e'ejuatie)ns  take'  the  fe)rm 
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in  a  ste'adily  rotating  frame'  ’..here'  71;  is  the  me'an  ve'leicity,  p  is  the'  meeelilie'd  me'an  |)re'ssure'. 
f?,  is  the  angular  velejcity  ejf  the  refe're'iiee'  frame',  r,,  =  i/'i/'  is  t  he'  Re'ynohls  si  re'ss  le'iiseer. 
anel  ;/  is  the  kine'inatic  vise-e)sity  e)f  the’  fluiel.  in  (l)-(2).  t  he’  Kinste’in  summatieen  ceenve'iit  ieeii 
applie's  te)  re’i)eate’d  inelice’s  and  e,,;.  eU’neete’s  the'  |)e’rmul at ieeti  le'nseer. 

First,  we  will  ce)nsiele’r  the  case’  eef  a  straight  e  ire  ular  pijee’  eer  neeii-cire  ular  duel  whe>s<’  axis 
lie’s  ale)ng  Mu'  :  diif'f’Me)n  ane!  whejse’  cre)ss-se’e  t ieui  lie's  in  the’  i/-piane'  (se’e’  Figure’  I).  Feir 


fully-developed  turbulent  flow,  where  u  =  u{x,y)i  +  u(.r,y)j  -f  w(.r,  y)k.  the  secondary  flow 
(u,  F)  is  derivable  from  a  stream  function  V’- 


The  Reynolds-averaged  Navier-Stokes  equations  can  then  be  solved  in  the  alternative  axial 
velocity/vorticity-streaTii  function  form  given  by 


_dw  _dw 
n— — 1-  V-—  —  G  + 
ax  ay 
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^  +  2n~ 

dx^  ^  dy 


V'V  =  C  (6) 

where  dpjdz  =  —G  is  the  constant  applied  pressure  gradient  driving  the  flow  and  = 
dvjdx—duldy  is  the  axial  component  of  the  mean  vorticity  vector.  Here,  the  angular  velocity 
fl  =  flj  corresponds  to  a  general  spanwise  rotation  since  we  can  align  the  component  of 
the  axis  of  rotation,  that  is  normal  to  the  axial  direction,  with  the  j/-axis.  Since  the  flow  is 
fully-developed,  the  component  of  the  angular  velocity  along  the  ^-axis  does  not  enter  into 
Eqs.  (4)-(5). 

From  (3)  and  (6),  it  is  clear  that  secondary  flows  are  generated  by  the  axial  mean  vorticity 
(  which  becomes  non-zero  in  a  straight  duct  or  pipe  only  if  there  are  normal  Reynolds  stress 
differences  {vyy  —  t^x)  or  Coriolis  forces  {2ildw f dy)  arising  from  a  spanwise  rotation.  In  the 
subsections  to  follow  we  will  briefly  categorize  secondary  flows  in  pipes  and  ducts. 

2.1.  Straight  Ducts  (17  =  o) 

It  is  the  axial  vorticity  source  term 

^  O^i'^TX  ~~  Tyy)  ^^Gy  ^^'^xy 

=  dtdy  + 

that  leads  to  the  generation  of  turbulent  secondary  flows  in  straight  non-circular  ducts.  A 
necessary  condition  for  the  occurrence  of  secondary  flows  is  that  the  axial  mean  velocity 
w{x,y)  must  give  ri.se  to  a  non-zero  normal  Reynolds  stress  difference  Txx  —  Tyy  (otherwise, 
Txy,  and  hence  will  vanish;  .see  Speziale  (I2,)()J).  For  a  circular  pipe,  even  though  —  Tyy 
is  non-zero,  <!>(-=  0  due  to  the  azimuthal  symmetry  of  the  axial  mean  velocity;  consecpient  ly 
no  secondary  flows  occur  [3].  For  non-circular  ducts,  if  —  Tyy  ^  0,  then  invariably,  4)(;  will 


be  non-zero  and  secondary  flows  are  generated.  From  tliis  result  it  is  rlt'ar  that  any  ('(idy 
viscosity  model  based  on  the  Boussinesq  hypothesis 


TT—  +  TT- 


(where  K  =  |r„  is  the  turbulent  kinetic  energy  and  //j  is  the  eddy  viscosity)  will  be  incapable 
of  predicting  secondary  flows  in  a  straight  non-circular  duct.  According  to  (8),  u  =  (O.O.Tr) 
yields  Tyy  =  t^x  —  |A'  and,  hence,  the  axial  mean  velocity  gives  rise  to  a  vanishing  normal 
Reynolds  stress  difference  Tyy  —  r^x  which  violates  the  necessary  condition  for  secondary 
flow.  Hence,  anisotropic  eddy  viscosity  models  -  where  yionlinear  strain-dependent  terms  are 
included  -  constitute  the  simplest  level  of  Reynolds  stress  closure  that  can  predict  secondary 
flows  in  straight  non-circular  ducts.  Three  examples  of  anisotropic  eddy  viscosity  models  are 
the  nonlinear  K  —  £  model  of  Speziale  [22],  the  two-scale  DIA  model  of  Yoshizawa  [23],  and 
the  RNG  based  model  of  Rubinstein  and  Barton  [24].  In  the  former  model  -  for  which  samjjle 
computations  will  be  presented  in  the  next  section  -  the  Reynolds  stress  is  represented  as 
follows: 

Aj  =  ~  ~  kehij\  (9) 


where 


dSi 


ot  dxk 


_ 

dxk  ^ 


is  the  frame-indifferent  Oldroyd  derivativeof  the  mean  rate  of  strain  tensor  Sij  =  ~[duildxj-\- 
duj/dxi)  and  e  is  the  turbulent  dissipation  rate;  =  0.09  and  Co  =  1-68  are  empirical 
constants.  In  the  limit  as  Co  0,  the  eddy  viscosity  relation  of  the  standard  K  —  £  model  is 
recovered.  It  can  be  shown  that  in  fully-developed  duct  flow,  the  axial  mean  velocity  w  gives 
rise  to  the  non-zero  normal  Reynolds  stress  difference  Tyy  —  Txx  =  C  oC'ffK^  j  £^)[{div  j  d.rY  — 
[dw/dyY]  in  (9).  Hence,  secondary  flows  are  generated  by  the  nonlinear  K  —  £  model. 

For  developing  secondary  flows,  where  history  effects  are  important,  a  full  Reynolds  stress 
closure  is  needed  for  a  more  complete  description.  Second-order  closures  are  based  on  the 
Reynolds  stress  transport  equation  [25]: 


dT,j  ,  ^  dTi, 


duj  dui  2 

dxk  dxk  3 


4“  //V'  Tij  2fl„t (c -|-  e.,„^iTj^^ 

where  H^j  is  the  combination  of  the  pressure-strain  correlation  and  the  deviatoric  part  of 
the  dissipation  rate  tensor  whereas  Vjj  is  the  turbulent  transport  term  (here,  1?  =  0  for  a 
stationary  duct).  From  (11),  it  is  straightforward  to  show  ~  by  substituting  u  =  (0,0.Tc) 


into  this  equation  -  that  secondary  flows  are  generated  by  0,^  and  which  are  the  only 
terms  that  yield  a  normal  Reynolds  stress  difference  Tyy  —  T:ci  (see  Speziale  [12,10]).  An  ordc-r 
of  magnitude  analysis,  at  high  Reynolds  numbers,  tends  to  indicate  that  anisotropi(>s  in  the 
pressure-strain  correlation  are  predominantly  responsible  for  the  generation  of  secondary 
flows;  anisotropies  in  the  turbulent  dissipation  rate  and  turbulent  transport  term  appear  to 
play  a  smaller  role. 


2.2.  Rotating  Pipes  and  Ducts 


In  a  rotating  circular  pipe,  the  axial  mean  velocity  uJ(r)  gives  rise  to  a  vanishing  'J>(;  by 
symmetry  arguments;  hence,  secondary  flows  are  generated  by  the  C’oriolis  term  'IQdTUfdy 
alone  (the  normal  Reynolds  stress  differences  in  only  have  an  indirect  effect  in  determining 
the  structure  of  the  resulting  fully-developed  secondary  flow).  Since  the  secondary  flows  in 
rotating  circular  pipes  are  generated  exclusively  by  Coriolis  forces,  eddy  viscosity  models  such 
as  the  K  —  e  model  are  capable  of  describing  this  effect  [9],  albeit  after  some  modifications. 
For  more  detailed  descriptions  of  the  flow  -  or  for  the  developing  flow  case  -  either  anisotropic 
eddy  viscosity  models  or  second-order  closures  should  be  used. 

Secondary  flows  in  rotating  non-circular  ducts  are  generated  by  two  sources:  normal 
Reynolds  stress  differences  embodied  in  the  term  and  Coriolis  forces  represented  by  the 
term  2Qdwldy.  Hence,  the  simplest  models  that  will  yield  acceptable  predictions  of  this 
flow  are  two-equation  models  with  an  anisotropic  eddy  viscosity.  However,  since  the  Coriolis 
forces  have  a  direct  effect  on  the  evolution  of  the  Reynolds  stresses  (see  Eq.  (11)),  second- 
order  closure  models  are  needed  for  a  more  complete  description  of  this  flow.  Even  in  the 
absence  of  secondary  flows,  the  Coriolis  forces  in  (11)  cause  the  axial  mean  velocity  profiles 
to  become  asymmetric  -  an  effect  that  is  difficult  to  describe  with  two-equation  models. 


2.3.  Curved  Pipes  and  Ducts 


Secondary  flows  in  curved  pipes  and  ducts  are  generated  by  centrifugal  forces.  This  can 
be  easily  seen  for  the  case  of  fully-developed  curved  duct  flow  with  the  streamwise  mean 
velocity  w  and  secondary  flow  (corresponding  to  the  directions  0,i\,  and  c,  respectively). 
Here,  the  mean  flow  equations  can  be  written  in  the  form: 


u 


dw 

dr 


u  dw  u  w  xu\  Id,  , 

-  +  T—  -f -  =  c;  d-  //  ^  -  -TrirTrff) 

oz  V  \  /  r  or 


dr.g 

dz 


(12) 

V 


where  (  =  (VXu)^  =  Oufdz  —  dvidr  is  the  streaniwise  mean  vortirity,  (/’  is  the  sc'ceiidary 
flow  stream  function,  and  —G  is  the  applied  mean  pressun*  gradient.  It  is  clear  from  (1.1) 
that  the  streaniwise  vorticity  source  term  {'2Wjr)dwldz  ~  which  constitutes  a  cc'ntrifugal 
acceleration  term  generated  by  the  primary  flow  tiJ-  is  the  main  generator  of  si'condary  flows 
in  curved  pipes.  In  curved  non-circular  ducts,  secondary  flows  are  generated  by  centrifugal 
effects  as  well  as  by  the  normal  Reynolds  stress  differences  —  T^r  and  T00  —  r^r  which  play  a 
crucial  role  when  the  curvature  is  weak,  (fon.sequently,  two-equation  models  with  an  isotropic 
eddy  viscosity,  such  as  the  K  —  e  model,  have  yielded  reasonably  acceptable  j)redictions  for 
fully-developed  secondary  flows  in  curved  pipes  and  curved  d>icts  with  modc'iate  to  strong 
curvature  ratios  [17].  In  order  to  analyze  curved  duct  flows  for  a  range  of  curvature  ratios, 
or  under  developing  conditions,  anisotropic  eddy  viscosity  models  or  second-order  closure 
models,  respectively,  are  needed. 

3.  ILLUSTRATIVE  EXAMPLES 

We  will  first  present  computations  of  fully-developed  turbulent  flow  in  a  straight  rectan¬ 
gular  duct  (see  Figure  1),  A  secondary  flow  with  eight  counter-rotating  vortices  is  generated 
by  normal  Reynolds  stress  differences.  Since  the  flow  is  fully-developed,  a  two-c'quation 
model  with  an  anisotropic  eddy  viscosity  -  solved  in  conjunction  with  wall  functions  suf¬ 
fices.  In  Figure  2,  the  mean  secondary  flow  streamlines  predicted  by  the  nonlinear  K  —  t 
model  (9)  at  high  Reynolds  numbers  are  shown  for  a  square  duct.  This  secondary  flow  is  of 
the  order  of  1-2%  of  the  axial  mean  velocity.  The  characteristic  eight- vortex  secondary  flow 
illustrated  in  the  schematic  provided  in  Figure  1  is  reproduced.  These  results  stand  in  sharp 
contrast  to  those  obtained  from  the  standard  K  —  e  model  which  erroneously  predicts  a  mu- 
directional  mean  turbulent  flow  -  a  deficiency  that  arises  from  the  use  of  an  eddy  viscosity 
model  based  on  the  Boussinesq  hypothesis  (8). 

More  quantitative  comparisons  will  now  be  made  for  a  straight  8  x  1  rectangular  duct.  The 
secondary  flow  streamlines  and  contours  of  the  normal  Reynolds  stress  anisotropy  Tj-j.  — 


predicted  by  the  nonlinear  K  —  e  model  are  compared  with  the  experimental  data  of  lloagland 
[1]  in  Figures  3(a)-3(b).  These  results  are  of  a  comparable  level  of  quality  as  those  obtaiiu'd 
from  a  second-order  closure  model  similar  to  that  of  Launder,  Reece  and  Rodi  [25]  as  shown  in 
Figures  4(a)-4(b).  Hence,  it  is  clear  that  two-equation  turbulence  moth-ls  with  an  anisotrojjic 
eddy  viscosity  yield  acceptable  predictions  for  fully-develoju'd  secondary  flows  in  a  straight 
rectangular  duct.  For  developing  turbulent  flows  in  non-circular  ducts,  s«-coiid-ordc'r  closun' 
models  yield  more  complete  predictions  [26]. 

The  fully-developed  secondary  flow  in  a  curve<l  s(|uare  duct  predicted  by  the  nonlin¬ 
ear  I\  —  £  model  at  high  Reynolds  numbers  is  shown  in  Figures  .5(a)-.5(b)  for  two  different 
curvature  ratios  Cr  (i-f'-,  the  ratio  of  the  radius  of  curvature  to  the  duct  width).  For  mod¬ 
erate  curvature,  there  is  a  double-vortex  secondary  flow  that  undergoes  a  liifurcation  tt)  a 
four-vortex  secondary  flow  when  the  curvature  ratio  TV  ~  10.  This  is  analogous  to  the 
Clortler  instability  in  curved  channel  flow  [27],  For  extremely  weak  curvature  >  10'^). 
the  double-vortex  secondary  flow  generated  by  centrifugal  effects  interacts  with  the  eight- 
vortex  secondary  flow  generated  by  normal  Reynolds  stress  differences  yielding  an  extremely 
complex  secondary  flow  pattern  which  we  will  not  show  herein  (.see  Hur  et  al.  [20]  and  Hur 
[28]  for  more  details).  An  anisotropic  eddy  viscosity  model  is  the  simplest  type  of  model 
that  can  describe  the  full  range  of  curvature  ratios.  For  a  more  complete  description  of  the 
flow  -  especially  under  developing  conditions  -  second-order  closure  models  are  preferable. 
However,  it  is  interesting  to  note  that  reasonably  acceptable  mean  flow  predictions  have  been 
obtained  for  fully-developed  curved  pipe  flow  using  the  standard  A’  —  t  model  for  a  limited 
range  of  curvature  ratios  [5]. 

Now,  we  will  discuss  the  prediction  of  curved  turbulent  pi])e  flows  that  are  not  fully- 
developed,  namely,  the  case  of  a  circular  pipe  f^-bend.  Here,  both  history  and  near- wall 
effects  play  a  role;  con.sequently,  a  .second-order  closure  with  a  lu'ar-wall  turbulence  model 
constitutes  the  preferred  approach.  We  will  show  some  illustrative  computations  for  the 
Launder,  Reece  and  Rodi  model  with  the  near  wall  turbulence  model  of  Lai  and  So  [29]. 
In  Figures  6(a)-6(b),  the  computed  .secondary  flow  and  mean  velocity  profiles  along  the 
pipe  centerline  (AA)  and  vertical  radius  (BB)  are  shown  to  compare  favorably  with  the 
experimental  data  of  Anwer  et  al.  [19]  for  a  location  67.5°  into  the  f^-bend.  At  this  same 
location  {9  —  67.5°),  the  secondary  flow  field  predicted  by  the  full  second-order  closure  model 
is  compared  with  its  counterpart  obtained  froiTi  the  A'  —e  model  in  Figures  7(a)-7(b).  It  is 
clear  from  these  results  that  the  second-order  closure  yields  a  more  detailed  picture  of  the 
secondary  flow  structure.  Ihdike  the  I\  —e  model,  the  second-order  closure  model  is  abh*  to 
])redict  the  i)resence  of  a  small  subsidiary  secoiulary  flow  cell  near  the  outer  bend  of  the  pipe 
(see  the  lower  h'ft-hand  corner  of  Figure  7(a)).  The  exist('nc('  of  this  sc'rondarv  flow  C('ll  has 


(i 


been  extrapolated  from  experiments  [21]. 

Finally,  we  will  present  some  eompute<l  rt'sidts  for  fiilly-devcdoped  t  nrbiilcMil  How  in  rei  l- 
angiilar  ducts  subjected  to  a  spanwise  rotation  (i.»'.,  t.lu'  duct  contig>irat ion  shown  in  I'igure 
1  mounted  in  a  frame  that  is  rotating  steadily  about  tlu'  //-axis  wit  h  an  angular  vcdocily 
n).  In  a  low-aspect-ratio  duct,  for  weak  to  moderat<'  rotations,  secondary  Hows  occur  that 
are  qualitatively  very  similar  to  those  obtained  in  curvc'd  rectangidar  ducts  (sec'  I’ignres 
5(a)-5(b)).  For  the  sake  of  brevity,  we  will  not  show  llu-sc'  double'- vortc'x  and  fonr-vortc'x 
secondary  flow  .solutions  (see  Younis  [HO]).  For  a  largc'-aspect-ral  io  rc'ct angular  duct  which 
is  used  to  experimentally  simulate  channel  flow  a  roll  instability  can  occur  at  intc'rmc'diatc' 
rotation  rates.  This  is  illustrated  in  Figure  8  which  shews  the  apix'arancc'  of  count c'r-rot at ing 
Taylor  vortices  in  the  interior  of  the'  duct  computc'd  using  a  nonliiu'ar  algc'braic  Kc'ynolds 
stress  model  for  a.  rotation  number  (i.e.,  angular  vc'locitv  normalizc'd  by  the'  channc'l  thicknc'ss 
and  bulk  mean  velocity  )  Ru  a;  0.1.  These  d'avlor  cc'lls  will  cause  the'  axial  mean  vc'locit  v 
along  the  duct  centerline  (which  is  n.se'd  to  approximate'  channc'l  How)  to  bc'coiTie'  (isiinimi  trie. 
Such  an  asymmetry  also  arises  in  the  absence'  of  se'condary  Hows  due'  to  t  he  dire'ct  e'lfc'ct  of 
CJoriolis  forces  on  the  Ke*ynolds  strc'ssc's  as  given  in  Fep  (11). 

One-dimensional  mean  ve'locity  calculations  of  rotating  channel  How  oblaine'd  using  the' 
SSG  second-order  closure  model  [fll]  with  wall  functions  are  compared  with  e'xpe'iime'iits 
[32]  in  Figures  9(a)-9(b).  At  we'ak  rotation  rates  (se'e  Figure'  9(a))  the  se'c'ond-order  closure' 
model  does  a  reasonably  good  job  of  pre'dictiiig  the'  a.symme'try  in  the  mc'an  ve'locity  profile'. 
However,  as  the  rotation  numl)er  Ro  becomes  of  the  order  of  10“' .  the  ejiiantit alive  ae  curacy 
of  the  residts  degrade  (se’e  Figure  9(b)  and  also  the  re'sidts  of  baemeh'r  e't  al.  [33]).  This 
could  be  partially  due  to  the  negle'ct  'f  roll  instabilitie's  which  have  be'en  doe  ume'iite'd  e'xpe'r- 
imentally  [8].  Full  second-order  closure  mode'l  calculations  of  a  rotating  re'ctangidar  duct  of 
large-aspect  ratio  -  as  illustrate'd  in  Figure  8  shoidel  be  conelucte'd  to  re'solve  this  issue'.  No 
such  detailed  calculations  of  roll-instabilitie's  in  rotating  turbulent  channel  flow  have'  ye't  to 
be  undertaken  with  a  second-order  closure. 

4.  CONCLUDING  REMARKS 

A  broad  overview  of  turbulent  se'cejiidarv  fle)ws  in  pipe's  anel  elucts  has  be'e'ii  pre'se'iite'd 
which  highlights  the  underlying  physical  nu'chanisms  re'sponsible  fe)r  their  generatiejii  anel  the' 
predictive  capabilities  of  Reynolds  stress  me)e|els  in  eh'seribing  tlu'se  fle)ws.  Se'conelarv  fle)ws 
arising  from  normal  Reynolds  stress  diffe’rene-es,  curvature,  anel  a  syste'in  re)tat  ie)n  have'  be'e'ii 
considered  in  an  effort  to  e'stablish  a  sufficie'iitly  ge'neral  basis  for  the  e'valuatie)n  of  me)ele'ls.  In 
the  opinion  of  the  authors,  two-equation  me>elels  with  an  anise)tre)pie'  e'dely  vise  e)sity  re'pre'se'iit 
the  simplest  level  of  model  that  can  pre'elict  a  wide  range  of  the'se  fle)ws  withe)ut  the'  ad  he)c 
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adjustment  of  constants  or  the  ad  hoc  prescription  of  turbulent  length  and  time  scales.  These 
two-equation  models  do  a  reasonably  good  job  of  predicting  turbulent  sc'condary  flows  in  their 
fully-developed  state.  For  a  mor»"  complete  description  of  these  flows  particularly  under 
developing  conditions  where  history  effects  or  body  forces  play  a  significant  role  st’cond- 
order  closures,  with  an  asymptotically  consistent  near-wall  turbulence  model,  are  j)ref('rred. 
While  future  research  is  still  needed,  the  results  presented  in  this  paper  demonstrate  the 
considerable  progress  that  has  been  made  during  the  past  two  decades  in  the  |)rediction  of 
turbulent  secondary  flows. 
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Figure  3.  Comparison  of  the  predictions  of  the  nonlinear  l\  —z  model  with  the  experimental 
data  of  Hoagland  [1]  for  3  x  1  rectangular  duct:  (a)  secondary  flow  streamlines  and  (h) 
contours  of  the  normal  Reynolds  stress  difference  —  r,^y. 
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Figure  4.  Comparison  of  the  predictions  of  a  second-order  closure  model  with  the  experi¬ 
mental  data  of  Hoagland  [1]  for  a  3  x  1  rectangular  duct:  (a)  secondary  flow  streamlines  and 
(b)  contours  of  the  normal  Reynolds  stress  difference  t^x  —  Tyy. 
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Figure  6.  Comparison  of  the  mean  velocity  predictions  of  the  second-order  closure  model 
of  Lai  and  So  [29]  with  experimental  data  [19]  for  a  circular  pipe  t/-bend  {0  =  67.5°):  (a) 
profiles  along  the  pipe  centerline  (AA)  and  (b)  profiles  along  the  vertical  radius  (BB). 
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Figure  9.  Comparison  of  the  mean  velocity  profiles  in  rotating  channel  flow.  (  SSCJ 
second-order  closure  model;  □  Experimental  data  {32]):  (a)  Ho  =  0.068.  Rx  =  .3.'), 000  and 
(b)  Ro  =  0.21,  Re  =  11,500. 
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